Genome-Wide Analysis of the Polygalacturonase Gene Family Sheds Light on the Characteristics, Evolutionary History, and Putative Function of Akebia trifoliata.

Polygalacturonase (PG) is one of the largest families of hydrolytic enzymes in plants. It is involved in the breakdown of pectin in the plant cell wall and even contributes to peel cracks. Here, we characterize PGs and outline their expression profiles using the available reference genome and transcriptome of Akebia trifoliata. The average length and exon number of the 47 identified AktPGs, unevenly assigned on 14 chromosomes and two unassembled contigs, were 5399 bp and 7, respectively. The phylogenetic tree of 191 PGs, including 47, 57, 51, and 36 from A. trifoliata, Durio zibethinus, Actinidia chinensis, and Vitis vinifera, respectively, showed that AktPGs were distributed in all groups except group G and that 10 AktPGs in group E were older, while the remaining 37 AktPGs were younger. Evolutionarily, all AktPGs generally experienced whole-genome duplication (WGD)/segmental repeats and purifying selection. Additionally, the origin of conserved domain III was possibly associated with a histidine residue (H) substitute in motif 8. The results of both the phylogenetic tree and expression profiling indicated that five AktPGs, especially AktPG25, could be associated with the cracking process. Detailed information and data on the PG family are beneficial for further study of the postharvest biology of A. trifoliata.

Exploring the Profile Contributions in Meyerozyma guilliermondii YB4 under Different NaCl Concentrations Using GC-MS Combined with GC-IMS and an Electronic Nose.

Using Meyerozyma guilliermondii YB4, which was isolated and screened from southern Sichuan pickles in the laboratory, as the experimental group, we investigated the changes in growth, total ester content, and volatile flavor substances of M. guilliermondii YB4 under different NaCl concentrations. The growth of M. guilliermondii YB4 was found to be inhibited by NaCl, and the degree of inhibition increased at higher NaCl concentrations. Additionally, the total ester content of the control group (CK) was significantly lower compared to the other groups (p < 0.05). The application of NaCl also resulted in distinct changes in the volatile profile of YB4, as evidenced by E-nose results. Gas chromatography-mass spectrometry (GC-MS) and gas chromatography-ion mobility spectrometry (GC-IMS) were employed to analyze the volatile compounds. A total of 148 and 86 volatiles were detected and identified using GC-MS and GC-IMS, respectively. Differential volatiles among the various NaCl concentrations in YB4 were determined by a variable importance in projection (VIP) analysis in partial least squares-discriminant analysis (PLS-DA). These differentially expressed volatiles were further confirmed by their relative odor activity value (ROAV) and odor description. Ten key contributing volatiles were identified, including ethanol, 1-pentanol, nonanal, octanal, isoamyl acetate, palmitic acid ethyl ester, acrolein, ethyl isobutanoate, prop-1-ene-3,3'-thiobis, and 2-acetylpyrazine. This study provides insights into the specificities and contributions of volatiles in YB4 under different NaCl concentrations. These findings offer valuable information for the development of aroma-producing yeast agents and the subsequent enhancement in the flavor of southern Sichuan pickles.

Identification Strategy of Biological Half Sibling Relationship.

OBJECTIVES: To compare the application value of the likelihood ratio (LR) method and identity by state (IBS) method in the identification involving half sibling relationships, and to provide a reference for the setting of relevant standards for identification of half sibling relationship. METHODS: (1) Based on the same genetic marker combinations, the reliability of computer simulation method was verified by comparing the distributions of cumulated identity by state score (CIBS) and combined full sibling index in actual cases with the distributions in simulated cases. (2) In different numbers of three genetic marker combinations, the simulation of full sibling, half sibling and unrelated individual pairs, each 1 million pairs, was obtained; the CIBS, as well as the corresponding types of cumulative LR parameters, were calculated. (3) The application value of LR method was compared with that of IBS method, by comparing the best system efficiency provided by LR method and IBS method when genetic markers in different amounts and of different types and accuracy were applied to distinguish the above three relational individual pairs. (4) According to the existing simulation data, the minimum number of genetic markers required to distinguish half siblings from the other two relationships using different types of genetic markers was estimated by curve fitting. RESULTS: (1) After the rank sum test, under the premise that the real relationship and the genetic marker combination tested were the same, there was no significant difference between the simulation method and the results obtained in the actual case. (2) In most cases, under the same conditions, the system effectiveness obtained by LR method was greater than that by IBS method. (3) According to the existing data, the number of genetic markers required for full-half siblings and half sibling identification could be obtained by curve fitting when the system effectiveness reached 0.95 or 0.99. CONCLUSIONS: When distinguishing half sibling from full sibling pairs or unrelated pairs, it is recommended to give preference to the LR method, and estimate the required number of markers according to the identification types and the population data, to ensure the identification effect.

Evolutionary conservation of the fidelity of transcription.

Accurate transcription is required for the faithful expression of genetic information. However, relatively little is known about the molecular mechanisms that control the fidelity of transcription, or the conservation of these mechanisms across the tree of life. To address these issues, we measured the error rate of transcription in five organisms of increasing complexity and found that the error rate of RNA polymerase II ranges from 2.9 × 10-6 ± 1.9 × 10-7/bp in yeast to 4.0 × 10-6 ± 5.2 × 10-7/bp in worms, 5.69 × 10-6 ± 8.2 × 10-7/bp in flies, 4.9 × 10-6 ± 3.6 × 10-7/bp in mouse cells and 4.7 × 10-6 ± 9.9 × 10-8/bp in human cells. These error rates were modified by various factors including aging, mutagen treatment and gene modifications. For example, the deletion or modification of several related genes increased the error rate substantially in both yeast and human cells. This research highlights the evolutionary conservation of factors that control the fidelity of transcription. Additionally, these experiments provide a reasonable estimate of the error rate of transcription in human cells and identify disease alleles in a subunit of RNA polymerase II that display error-prone transcription. Finally, we provide evidence suggesting that the error rate and spectrum of transcription co-evolved with our genetic code.

The PPR-Domain Protein SOAR1 Regulates Salt Tolerance in Rice.

Previous studies in Arabidopsis reported that the PPR protein SOAR1 plays critical roles in plant response to salt stress. In this study, we reported that expression of the Arabidopsis SOAR1 (AtSOAR1) in rice significantly enhanced salt tolerance at seedling growth stage and promoted grain productivity under salt stress without affecting plant productivity under non-stressful conditions. The transgenic rice lines expressing AtSOAR1 exhibited increased ABA sensitivity in ABA-induced inhibition of seedling growth, and showed altered transcription and splicing of numerous genes associated with salt stress, which may explain salt tolerance of the transgenic plants. Further, we overexpressed the homologous gene of SOAR1 in rice, OsSOAR1, and showed that transgenic plants overexpressing OsSOAR1 enhanced salt tolerance at seedling growth stage. Five salt- and other abiotic stress-induced SOAR1-like PPRs were also identified. These data showed that the SOAR1-like PPR proteins are positively involved in plant response to salt stress and may be used for crop improvement in rice under salinity conditions through transgenic manipulation.

Multiomics analysis elucidated molecular mechanism of aromatic amino acid biosynthesis in Akebia trifoliata fruit.

Akebia trifoliata is a novel edible and healthy fruit. Here, we found that this fruit had the highest content of total free amino acids and three aromatic amino acids (AAAs) compared with the other popular fruits, and there was an obvious inverse relationship between AAA and flavonoid levels in various fruit tissues. Multiomics analysis revealed that the evolutionarily strengthened synthetic pathway of all three AAAs, the largely regulating ability conferred by ASP5 in the arogenate pathway and the complementary phenylpyruvate pathway endorsed by ADT of both Phe and Tyr biosynthesis provided reasonable explanations for the high AAA content in the flesh of A. trifoliata fruit. Gene-specific expression could be the main reason for the inverse relationship between AAAs and flavonoids. This study will help us understand the metabolic mechanism of AAAs and to develop A. trifoliata as a fresh fruit crop and medicinal plant by molecular breeding strategies.

The chromosome-level genome of Akebia trifoliata as an important resource to study plant evolution and environmental adaptation in the Cretaceous.

The environmental adaptation of eudicots is the most reasonable explanation for why they compose the largest clade of modern plants (>70% of angiosperms), which indicates that the basal eudicots would be valuable and helpful to study their survival and ability to thrive throughout evolutionary processes. Here, we detected two whole-genome duplication (WGD) events in the high-quality assembled Akebia trifoliata genome (652.73 Mb) with 24 138 protein-coding genes based on the evidence of intragenomic and intergenomic collinearity, synonymous substitution rate (KS ) values and polyploidization and diploidization traces; these events putatively occurred at 85.15 and 146.43 million years ago (Mya). The integrated analysis of 16 species consisting of eight basal and eight core eudicots further revealed that there was a putative ancient WGD at the early stage of eudicots (temporarily designated θ) at 142.72 Mya, similar to the older WGD of Akebia trifoliata, and a putative core eudicot-specific WGD (temporarily designated ω). Functional enrichment analysis of retained duplicate genes following the θ event is suggestive of adaptation to the extreme environment change in both the carbon dioxide concentration and desiccation around the Jurassic-Cretaceous boundary, while the retained duplicate genes following the ω event is suggestive of adaptation to the extreme droughts, possibly leading to the rapid spread of eudicots in the mid-Cretaceous. Collectively, the A. trifoliata genome experienced two WGD events, and the older event may have occurred at the early stage of eudicots, which likely increased plant environmental adaptability and helped them survive in ancient extreme environments.

Characterization of the MADS-Box Gene Family in Akebia trifoliata and Their Evolutionary Events in Angiosperms.

As the largest clade of modern plants, flower plants have evolved a wide variety of flowers and fruits. MADS-box genes play key roles in regulating plant morphogenesis, while basal eudicots have an evolutionarily important position of acting as an evolutionary bridge between basal angiosperms and core eudicots. Akebia trifoliata is an important member of the basal eudicot group. To study the early evolution of angiosperms, we identified and characterized the MADS-Box gene family on the whole-genome level of A. trifoliata. There were 47 MADS-box genes (13 type I and 34 type II genes) in the A. trifoliata genome; type I genes had a greater gene length and coefficient of variation and a smaller exon number than type II genes. A total of 27 (57.4%) experienced whole or segmental genome duplication and purifying selection. A transcriptome analysis suggested that three and eight genes were involved in whole fruit and seed development, respectively. The diversification and phylogenetic analysis of 1479 type II MADS-box genes of 22 angiosperm species provided some clues indicating that a γ whole genome triplication event of eudicots possibility experienced a two-step process. These results are valuable for improving A. trifoliata fruit traits and theoretically elucidating evolutionary processes of angiosperms, especially eudicots.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Akebiatrifoliata: A Bioinformatics Study.

WRKY transcription factors have been found in most plants and play an important role in regulating organ growth and disease response. Outlining the profile of WRKY genes is a very useful project for studying morphogenesis and resistance formation. In the present study, a total of 63 WRKY genes consisting of 13 class I, 41 class II, and 9 class III genes were identified from the newly published A. trifoliata genome, of which 62 were physically distributed on all 16 chromosomes. Structurally, two AkWRKY genes (AkWRKY6 and AkWRKY52) contained four domains, and AkWRKY17 lacked the typical heptapeptide structure. Evolutionarily, 42, 16, and 5 AkWRKY genes experienced whole genome duplication (WGD) or fragmentation, dispersed duplication, and tandem duplication, respectively; 28 Ka/Ks values of 30 pairs of homologous genes were far lower than 1, while those of orthologous gene pairs between AkWRKY41 and AkWRKY52 reached up to 2.07. Transcriptome analysis showed that many of the genes were generally expressed at a low level in 12 fruit samples consisting of three tissues, including rind, flesh, and seeds, at four developmental stages, and interaction analysis between AkWRKY and AkNBS genes containing W-boxes suggested that AkWRKY24 could play a role in plant disease resistance by positively regulating AkNBS18. In summary, the WRKY gene family of A. trifoliata was systemically characterized for the first time, and the data and information obtained regarding AkWRKY could be very useful in further theoretically elucidating the molecular mechanisms of plant development and response to pathogens and practically improving favorable traits such as disease resistance.

Identification and Characterization of NBS Resistance Genes in Akebia trifoliata.

Akebia trifoliata is an important multiuse perennial plant that often suffers attacks from various pathogens due to its long growth cycle, seriously affecting its commercial value. The absence of research on the resistance (R) genes of A. trifoliata has greatly limited progress in the breeding of resistant varieties. Genes encoding proteins containing nucleotide binding sites (NBSs) and C-terminal leucine-rich repeats (LRRs), the largest family of plant resistance (R) genes, are vital for plant disease resistance. A comprehensive genome-wide analysis showed that there were only 73 NBS genes in the A. trifoliata genome, including three main subfamilies (50 coiled coil (CC)-NBS-LRR (CNL), 19 Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) and four resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) genes). Additionally, 64 mapped NBS candidates were unevenly distributed on 14 chromosomes, most of which were assigned to the chromosome ends; 41 of these genes were located in clusters, and the remaining 23 genes were singletons. Both the CNLs and TNLs were further divided into four subgroups, and the CNLs had fewer exons than the TNLs. Structurally, all eight previously reported conserved motifs were identified in the NBS domains, and both their order and their amino acid sequences exhibited high conservation. Evolutionarily, tandem and dispersed duplications were shown to be the two main forces responsible for NBS expansion, producing 33 and 29 genes, respectively. A transcriptome analysis of three fruit tissues at four developmental stages showed that NBS genes were generally expressed at low levels, while a few of these genes showed relatively high expression during later development in rind tissues. Overall, this research is the first to identify and characterize A. trifoliata NBS genes and is valuable for both the development of new resistant cultivars and the study of molecular mechanisms of resistance.

BGISEQ-500RS sequencing of a 448-plex SNP panel for forensic individual identification and kinship analysis.

Next generation sequencing (NGS)-based single nucleotide polymorphism (SNP) genotyping is widely used in the field of forensics. SNP genotyping data from several NGS platforms have been published, but forensic application trials of DNA nanoball sequencing platforms have been very limited. In this work, we developed a 448-plex SNP panel on the BGISEQ-500RS platform. The sequencing metrics of a total of 261 samples that were sequenced with this panel are reported in detail. The average sequencing depth was 8373 × and the average heterozygosity of the 448-plex assay was 0.85. Sensitivity analysis showed that 325 SNPs were successfully genotyped with as little as 50 pg of genomic DNA, with the mean quality score of the sequencing data above Q30. Forensic parameters were calculated based on the data of 142 unrelated Chinese Han individuals and the combined matching probability was as low as 5.21 × 10-101. Kinship analyses based on experiments and computer simulations showed that the 448-panel was as effective as the ForenSeq™ DNA Signature Prep Kit for second-degree kinship identification, and when the two panels were merged, the related pairs were almost completely distinguished from unrelated pairs. The 448-plex SNP panel on the BGISEQ-500RS platform provides a powerful tool for forensic individual identification and kinship analysis.

Transcriptomic Analyses Shed Light on Critical Genes Associated with Bibenzyl Biosynthesis in Dendrobium officinale.

The Dendrobium plants (members of the Orchidaceae family) are used as traditional Chinese medicinal herbs. Bibenzyl, one of the active compounds in Dendrobium officinale, occurs in low amounts among different tissues. However, market demands require a higher content of thes compounds to meet the threshold for drug production. There is, therefore, an immediate need to dissect the physiological and molecular mechanisms underlying how bibenzyl compounds are biosynthesized in D. officinale tissues. In this study, the accumulation of erianin and gigantol in tissues were studied as representative compounds of bibenzyl. Exogenous application of Methyl-Jasmonate (MeJA) promotes the biosynthesis of bibenzyl compounds; therefore, transcriptomic analyses were conducted between D. officinale-treated root tissues and a control. Our results show that the root tissues contained the highest content of bibenzyl (erianin and gigantol). We identified 1342 differentially expressed genes (DEGs) with 912 up-regulated and 430 down-regulated genes in our transcriptome dataset. Most of the identified DEGs are functionally involved in the JA signaling pathway and the biosynthesis of secondary metabolites. We also identified two candidate cytochrome P450 genes and nine other enzymatic genes functionally involved in bibenzyl biosynthesis. Our study provides insights on the identification of critical genes associated with bibenzyl biosynthesis and accumulation in Dendrobium plants, paving the way for future research on dissecting the physiological and molecular mechanisms of bibenzyl synthesis in plants as well as guide genetic engineering for the improvement of Dendrobium varieties through increasing bibenzyl content for drug production and industrialization.

Insight-HXMT observations of jet-like corona in a black hole X-ray binary MAXI J1820+070.

A black hole X-ray binary produces hard X-ray radiation from its corona and disk when the accreting matter heats up. During an outburst, the disk and corona co-evolves with each other. However, such an evolution is still unclear in both its geometry and dynamics. Here we report the unusual decrease of the reflection fraction in MAXI J1820+070, which is the ratio of the coronal intensity illuminating the disk to the coronal intensity reaching the observer, as the corona is observed to contrast during the decay phase. We postulate a jet-like corona model, in which the corona can be understood as a standing shock where the material flowing through. In this dynamical scenario, the decrease of the reflection fraction is a signature of the corona's bulk velocity. Our findings suggest that as the corona is observed to get closer to the black hole, the coronal material might be outflowing faster.

Light Activates the Translational Regulatory Kinase GCN2 via Reactive Oxygen Species Emanating from the Chloroplast.

Cytosolic mRNA translation is subject to global and mRNA-specific controls. Phosphorylation of the translation initiation factor eIF2α anchors a reversible regulatory switch that represses cytosolic translation globally. The stress-responsive GCN2 kinase is the only known kinase for eIF2α serine 56 in Arabidopsis (Arabidopsis thaliana). Here, we show that conditions that generate reactive oxygen species (ROS) in the chloroplast, including dark-light transitions, high light, and the herbicide methyl viologen, rapidly activated GCN2 kinase, whereas mitochondrial and endoplasmic reticulum stress did not. GCN2 activation was light dependent and mitigated by photosynthesis inhibitors and ROS quenchers. Accordingly, the seedling growth of multiple Arabidopsis gcn2 mutants was retarded under excess light conditions, implicating the GCN2-eIF2α pathway in responses to light and associated ROS. Once activated, GCN2 kinase preferentially suppressed the ribosome loading of mRNAs for functions such as mitochondrial ATP synthesis, the chloroplast thylakoids, vesicle trafficking, and translation. The gcn2 mutant overaccumulated transcripts functionally related to abiotic stress, including oxidative stress, as well as innate immune responses. Accordingly, gcn2 displayed defects in immune priming by the fungal elicitor, chitin. Therefore, we provide evidence that reactive oxygen species produced by the photosynthetic apparatus help activate the highly conserved GCN2 kinase, leading to eIF2α phosphorylation and thus affecting the status of the cytosolic protein synthesis apparatus.

Joint and Direct Optimization for Dictionary Learning in Convolutional Sparse Representation.

Convolutional sparse coding (CSC) is a useful tool in many image and audio applications. Maximizing the performance of CSC requires that the dictionary used to store the features of signals can be learned from real data. The so-called convolutional dictionary learning (CDL) problem is formulated within a nonconvex, nonsmooth optimization framework. Most existing CDL solvers alternately update the coefficients and dictionary in an iterative manner. However, these approaches are prone to running redundant iterations, and their convergence properties are difficult to analyze. Moreover, most of those methods approximate the original nonconvex sparse inducing function using a convex regularizer to promote computational efficiency. This approach to approximation may result in nonsparse representations and, thereby, hinder the performance of the applications. In this paper, we deal with the nonconvex, nonsmooth constraints of the original CDL directly using the modified forward-backward splitting approach, in which the coefficients and dictionary are simultaneously updated in each iteration. We also propose a novel parameter adaption scheme to increase the speed of the algorithm used to obtain a usable dictionary and in so doing prove convergence. We also show that the proposed approach is applicable to parallel processing to reduce the computing time required by the algorithm to achieve convergence. The experimental results demonstrate that our method requires less time than the existing methods to achieve the convergence point while using a smaller final functional value. We also applied the dictionaries learned using the proposed and existing methods to an application involving signal separation. The dictionary learned using the proposed approach provides performance superior to that of comparable methods.

Adaptive ADMM for Dictionary Learning in Convolutional Sparse Representation.

In this paper, we propose a novel approach to convolutional sparse representation with the aim of resolving the dictionary learning problem. The proposed method, referred to as the adaptive alternating direction method of multipliers (AADMM), employs constraints comprising non-convex, non-smooth terms, such as the l0 -norm imposed on the coefficients and the unit-norm sphere imposed on the length of each dictionary element. The proposed scheme incorporates a novel parameter adaption scheme that enables ADMM to achieve convergence more quickly, as evidenced by numerical and theoretical analysis. In experiments involving image signal applications, the dictionaries learned using AADMM outperformed those learned using comparable dictionary learning methods.

Molecular and biochemical characterization of the jasmonic acid methyltransferase gene from black cottonwood (Populus trichocarpa).

Methyl jasmonate is a metabolite known to be produced by many plants and has roles in diverse biological processes. It is biosynthesized by the action of S-adenosyl-l-methionine:jasmonic acid carboxyl methyltransferase (JMT), which belongs to the SABATH family of methyltransferases. Herein is reported the isolation and biochemical characterization of a JMT gene from black cottonwood (Populus trichocarpa). The genome of P. trichocarpa contains 28 SABATH genes (PtSABATH1 to PtSABATH28). Recombinant PtSABATH3 expressed in Escherichia coli showed the highest level of activity with jasmonic acid (JA) among carboxylic acids tested. It was therefore renamed PtJMT1. PtJMT1 also displayed activity with benzoic acid (BA), with which the activity was about 22% of that with JA. PtSABATH2 and PtSABATH4 were most similar to PtJMT1 among all PtSABATHs. However, neither of them had activity with JA. The apparent Km values of PtJMT1 using JA and BA as substrate were 175µM and 341µM, respectively. Mutation of Ser-153 and Asn-361, two residues in the active site of PtJMT1, to Tyr and Ser respectively, led to higher specific activity with BA than with JA. Homology-based structural modeling indicated that substrate alignment, in which Asn-361 is involved, plays a role in determining the substrate specificity of PtJMT1. In the leaves of young seedlings of black cottonwood, the expression of PtJMT1 was induced by plant defense signal molecules methyl jasmonate and salicylic acid and a fungal elicitor alamethicin, suggesting that PtJMT1 may have a role in plant defense against biotic stresses. Phylogenetic analysis suggests that PtJMT1 shares a common ancestor with the Arabidopsis JMT, and functional divergence of these two apparent JMT orthologs has occurred since the split of poplar and Arabidopsis lineages.

Interlayer bit allocation for scalable video coding.

In this paper, we present a theoretical analysis of the distortion in multilayer coding structures. Specifically, we analyze the prediction structure used to achieve temporal, spatial, and quality scalability of scalable video coding (SVC) and show that the average peak signal-to-noise ratio (PSNR) of SVC is a weighted combination of the bit rates assigned to all the streams. Our analysis utilizes the end user's preference for certain resolutions. We also propose a rate-distortion (R-D) optimization algorithm and compare its performance with that of a state-of-the-art scalable bit allocation algorithm. The reported experiment results demonstrate that the R-D algorithm significantly outperforms the compared approach in terms of the average PSNR.

Developing consensus 3D-QSAR and pharmacophore models for several beta-secretase, farnesyl transferase and histone deacetylase inhibitors.

Three consensus 3D-QSAR (c-3D-QSAR) models were built for 38, 34, and 78 inhibitors of ß-secretase, histone deacetylase, and farnesyltransferase, respectively. To build an individual 3D-QSAR model, the structures of an inhibitor series are aligned through docking of a protein receptor into the active site using the program GOLD. CoMFA, CoMSIA, and Catalyst are then performed for the training set of each structurally aligned inhibitor series to obtain a 3D-QSAR model. Since the consensus in features identified is high for the same pharmacophore features selected for building a 3D-QSAR model by a 3D-QSAR method, a c-3D-QSAR model for each inhibitor series is constructed by combining the pharmacophore features selected for building the 3D-QSAR model using the SYBYL spread sheet and PLS module. Each c-3D-QSAR pharmacophore model built was examined visually and compared with that obtained by simultaneous mapping of the corresponding 3D-QSAR pharmacophores built onto a selected inhibitor structure. It was found that the c-3D-QSAR model built for an inhibitor series improves not only the overall prediction statistics for both training and test sets but also the prediction accuracy for some less active inhibitors of the series.

Biosynthesis and emission of insect-induced methyl salicylate and methyl benzoate from rice.

Two benzenoid esters, methyl salicylate (MeSA) and methyl benzoate (MeBA), were detected from insect-damaged rice plants. By correlating metabolite production with gene expression analysis, five candidate genes encoding putative carboxyl methyltransferases were identified. Enzymatic assays with Escherichia coli-expressed recombinant proteins demonstrated that only one of the five candidates, OsBSMT1, has salicylic acid (SA) methyltransferase (SAMT) and benzoic acid (BA) methyltransferase (BAMT) activities for producing MeSA and MeBA, respectively. Whereas OsBSMT1 is phylogenetically relatively distant from dicot SAMTs, the three-dimensional structure of OsBSMT1, which was determined using homology-based structural modeling, is highly similar to those of characterized SAMTs. Analyses of OsBSMT1 expression in wild-type rice plants under various stress conditions indicate that the jasmonic acid (JA) signaling pathway plays a critical role in regulating the production and emission of MeSA in rice. Further analysis using transgenic rice plants overexpressing NH1, a key component of the SA signaling pathway in rice, suggests that the SA signaling pathway also plays an important role in governing OsBSMT1 expression and emission of its products, probably through a crosstalk with the JA signaling pathway. The role of the volatile products of OsBSMT1, MeSA and MeBA, in rice defense against insect herbivory is discussed.